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ABSTRACT 

We report the polarization measurement in prompt 7-ray emission of GRB 100826A with the 
Gamma- Ray Burst Polarimeter (GAP) aboard the small solar power sail demonstrator IKAROS. We 
detected the firm change of polarization angle (PA) during the prompt emission with 99.9 % (3.5 a) 
confidence level, and the average polarization degree (n) of 27 ± 11 % with 99.4 % (2.9 a) confidence 
level. Here the quoted errors are given at 1 a confidence level for two parameters of interest. The sys- 
tematic errors have been carefully included in this analysis, unlike any previous reports. Such a high 
LT can be obtained in several emission models of gamma-ray bursts (GRBs), including synchrotron and 
photospheric models. However, it is difficult to explain the observed significant change of PA within 
the framework of axisymmetric jet as considered in many theoretical works. The non-axisymmctric 
(e.g., patchy) structures of the magnetic fields and/or brightness inside the relativistic jet are therefore 
required within the observable angular scale of ~ Our observation strongly indicates that the 

polarization measurement is a powerful tool to constrain the GRB production mechanism, and more 
theoretical works are needed to discuss the data in more details. 

Subject headings: gamma-ray: bursts — gamma-ray: polarization — gamma-ray: emission mechanism 



1. INTRODUCTION 

Gamma-ray bursts (GRBs) are the most energetic ex- 
plosions in the universe, and the isotropic luminosity 
reaches 10 54 erg s _1 for the brightest bursts. Since the 
discovery of X-ra y afterglow of GRBs by BeppoSAX 
(|Costa et al.lll997l) and the identification of optical coun- 
terparts at cosmological distance, many observational 
facts and theoretical works have led us to understand 
the nature of GRBs. However, a crucial issue still re- 
maining to answer is how to release such a huge energy 
as 7-ray photons. In spite of extensive discussions with 
the spectral and lightcurve information being collected, 
there are still several possible models of GRBs (Mcszaros 
2006). The polarimetric observations provide us com- 
pletely different information. A firm detection of linear 
polarization of 7 -ray photons will further constr ain the 
emission models (|Lazzatill200l iToma et al.ll2009| ). 

The polarization measurement of GRBs has been 
performed in the late optical afterglows at first 
(ICovino et al.l[l99l IWiiers et all Il999t IRol et alJ [2000l: 
iGorosabel et al.l 120041 : ICovino et al.l 120041 ). All of these 
result s show that the polar ization degree (n) of a few % 
level. iMundel et all (|2007[ ) set upper limit of n < 8 % 
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for the early optical afterglow. Detailed observations are 
performed for the optical afterglow of GRB 030329, and 
time variation s of the polarization degree and angle (PA) 
were reported ([Greiner et al.|[2CJ03l ). Recently, a polariza- 
tion degree of II = 10 ± 1 % was reported for the early 
optical afterglow of GR B 090102, at hund reds of seconds 
since the burst trigger ([Steele et alJl2009t ). 

On the other hand, there are controversies on polariza- 
tion detections in the prompt 7-ray emission of GRBs. 
The first report was the detection of a high linear po- 
larization of n = 80 ± 20 % at a c onfidence level of 
> 5.7 a by RHESSI from GRB 021206 (ICoburn fe Boggsl 
120031 ). although independent groups could not confirm 
any polarization signa l s by reanalyzing th e same data 
(|Rutledge fc Foxl 120041: iWigger et all I2004D . The sec- 
ond reports showed also high degr ee of polarization for 
GRB 041219A wit h n = 98 ±33 % (jKalemci et al.ll200l 
and n = 63±f % (jMcGlvnn et al.|[200l for the brightest 
pulse in the 100-350 keV energy band of INTEGRAL- 
SPI at 2 cr confidence level fo r 2 parameters of interest 
(n and PA). lGotz et~aTI ([20091 ) also reported more signifi- 
cant detection for time resolved analyses^, and suggested 
possible variable polarization properties for the same 
GRB 041219A in the 200-800 keV energy band obse rved 
with INTEGRAL-IBIS. However, iGotz et al.l ((20091) set 
a strict upper limit of n < 4 % for the same bright- 
est pulse showing high n in the INTEGR A L-SPI data 
(|Kalemci et alJl2007tlMcGlvnn et al.ll2007| ). ([Gotz et al.1 
(2009) explained that the null net polarization degree 
results from the superposition of signals with different 
polarization angles. However, this interpretation can- 
not explain the positive detection by INTEGRAL-SPI 
simultaneously.) These inconsistent results between SPI 

7 See Table 1 and Figure 3 of lGotz et all ||2009T ). One should note 
that their definition of confidence level is given for one parameter 
of interest. 
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and IBIS led to confusion again, and the existence of 
7-ray polarization is still in debate. As the authors 
themselves pointed out in their reports, these results are 
with low statistics of ~ 2 a level, and may be strongly 
affected by the instrumental systematics uncertainties 
(iKalemci et al.1 120071 : iMcGlvnn et~aT1 120071 iGotz et al.1 
2009). These controversies and conflicts indicate diffi- 
culties in detecting the 7-ray polarization. 

In this Letter, we report the detection of 7-ray polar- 
ization and also the change of PA for extremely bright 
GRB 100826A using the newly developed GRB polarime- 
ter. Our polarimcter is completely designed for the po- 
larization measurement of prompt GRBs, and well cal- 
ibrated during developing phase before launch. Espe- 
cially, using a detector of proto-flight model, we exper- 
imentally and numerically understand the response for 
polarized 7-ray with low systematic uncertainties. In the 
following sections, we show the observation (section 2), 
data analysis (section 3) and discussion for the emission 
mechanism of the prompt GRBs (section 4). 

2. OBSERVATIONS 

IKAROS (|Kawaguchi et al.1 |200H IMori et al.1 [2010^ 
(Interplanetary Kite-craft Accelerated by Radiation Of 
the Sun) is a small solar-power-sail demonstrator and 
was successfully launched on 21 May 2010. IKAROS has 
a large polyimide membrane of 20 m in diameter, and 
this translates the solar radiation pressure to the thrust 
of the spacecraft. Since the deployment of the sail on 
9 June 2010, IKAROS started the solar-sailing toward 
Venus. 

The Gamma-R ay B urst Polarimeter ( GAP) 
(lYonetoku et all 120111 : iMurakami et all l2010i 
I Yonetoku etaTT 120061) aboard IKAROS is fully de- 
signed to measure a degree of linear polarization in 
the prompt emission of GRBs in the energy range of 
70-300 keV. The detection principle is an anisotropy of 
Compton scattered photons. If the incident 7-rays are 
linearly polarized, the distribution of scattered photons 
is due to Klein-Nishina effect which approximately 
shows sin 2 <fi curves, where (f> is the scattering angle. 
The GAP consists of a central plastic scatterer of 17 cm 
in diameter and 6 cm in thickness, and the surrounding 
12 CsI(Tl) scintillators with 5 mm in thickness. The 
coincidence events within a gate time of 5 /its between 
the signal from any Csl and that from the plastic are 
accepted for the polarization analysis. The GAP's high 
axial symmetry in shape and the high gain uniformity 
are the keys for reliable measurement of polarizations, 
and to avoid a fake modulation due to background 
7-rays. There are no external parts of spacecraft inside 
the GAP field of view. Moreover the detector cases 
(chassis), except for the detector top, are covered by 
thin lead sheets with 0.5 mm thickness. Therefore the 
effect of background 7-rays scattered by the spacecraft 
body are negligible. 

The GAP detected GRB 100826A on 26 August 2010 
at 22:57:20.8 (UT) on the way to Venus at about 0.21 AU 
apart from the Earth. The lightcurve of the prompt emis- 
sion is shown in Figure [U This burst was also detected by 
several satellite s, and was localize d by an interplanetary 
network (IPN) (|Hurlev et alll2010l ). Combining the GAP 
data with the published IPN information, the direction of 
this burst is derived as (a, 5) = (279.6±0.2, -22.3±0.5), 



which corresponds to 20.0 degree off-axis from the cen- 
ter of the GAP field of view. An energy fluence of this 
burst is (3.0 ± 0.3) x 10~ 4 erg cm" 2 in 20 keV-10 MeV 
band by KONUS (|Golenetskii et all 120101 ) . which is the 
top 1 % of the brightest events listed in BATSE cata- 
log. The low- and high-energy photon indices are re- 
ported as olb — — 1.311q;q5 and (3b — ~2.1+q;J,, respec- 
tively, and the vF v peak energy as E p — 6061 109 keV 
(Golcnets kii et all [20101 An optical afterglow of this 
GRB was not reported, so its redshift is unknown. 

3. DATA ANALYSIS 

We divided the entire data into two time intervals as 
labeled Interval-1 and -2 in a lightcurve of Figured] The 
total numbers of 7-ray photons after subtracting back- 
ground are 32,924 and 19,007 photons for each interval, 
respectively. The first part of this burst shows a large 
flare lasting 47 seconds since the trigger, and a follow- 
ing 53 seconds consists of multiple spikes. Although the 
Interval-2 has several spikes, we combined all of them to 
keep photon statistics. 

We used the background modulation curve obtained 
by the 36.7 hours integration just before and after the 
GRB trigger time. In this case, the averaged coincidence 
background rate is quite stable as 5.6 counts s _1 CsI _1 , 
and its modulation curve can be described as constant 
with the standard deviation of 0.1 % level. After sub- 
tracting the background modulation, the total numbers 
of the coincidence 7-rays (polarization signals) are 4,821 
and 2,733 for Interval-1 and -2, respectively. 

3.1. Model Functions 

When we perform the polarization analysis, we cal- 
culated the model function (or the detector response) 
for the polarized gamma-rays with Geant 4 Monte-Carlo 
simulato r. A geometrical mass model of GAP was intro- 
duced in lYonetoku et al.1 (|2011[ ). We can set all the pos- 
sible situations with different off- axis angles, azimuthal 
phase angles, spectral parameters, as well as polarization 
degrees and angles of incoming 7-rays. In this analysis, 
the free parameters are the polarization degrees and an- 
gles, since the other parameters are constrained as shown 
in section 2. 

Simulating the interactions between the incoming 
gamma-rays and the geometrical mass model with 
Monte-Carlo method, we created the model functions of 
the modulation curves. Then we selected the coincidence 
events between the plastic and one Csl scintillator as the 
polarization signals. We simulated the model function 
with the step resolutions of 5 % for polarization degrees 
and 5 degrees for phase angles. We will fit the observed 
modulation curves with the model function by a least 
squares method in the following subsections. 

3.2. Systematic Uncertainty 

We should treat not only the statistical error but also 
the systematic uncertainties correctly. Here, the system- 
atic uncertainties mean any unexplainablc difference be- 
tween the observed modulation curve and the simulated 
one. 

Imperfect tunings of parameters in the ground and 
in-orbit calibrations for the incident radiation, espe- 
cially from off-axis direction, cause more important sys- 
tematic uncertainty. First, we performed some ground 
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based experiments with non-polarized radio isotopes 
57 Co (122 keV) and 241 Am (59.5 keV). They were set 
at the distance of 1.0 m from the center of proto- flight 
model with several incident angles from zero to 50 de- 
grees with the step resolutions of 5 degrees. We measured 
experimental modulation curves for each setup. After 
that, on the assumption of non-polarized 7-rays with the 
monochromatic energy of 122 keV and/or 59.5 keV, we 
calculated the numerical modulation curve with Geant 4 
simulator. 

Then, of course, the both modulation curves show 
some discrepancy which cannot be explained within the 
statistical uncertainty, because of instrumental system- 
atic uncertainties. We calculated the standard deviation 
of two modulations and derived the systematic errors to 
be <7 sys = 1.8 % of the total coincidence 7-rays for the 
case of 20.0 degree incident angle. Although this value 
is not negligible, the statistical error (about a sta t ~ 5 % 
leves) still dominates the systematic one in this case. We 
included the systematic uncertainty into the total errors 
as a totai = a stat + a sys f° r eacn bin of polarization data. 

3.3. Fitting the Polarization Data 

First of all, we performed the polarization analysis for 
the whole dataset of time interval of 0-100 seconds. Then 
we obtained an acceptable result with the model of non- 
polarized modulation curve, and set only upper limit 
of II < 30 % (2 a confidence level). We consider the 
possibility that the polarization degree may be weak if 
the polarization angle cha nges during t he en tire burst 
duration, as sugg ested by IGotz et all ([2009). There- 
fore, we performed time resolved polarization analyses 
for datasets of Interval- 1 and -2. In Figure [2j we show 
the observed source modulation curves as a function of 
scattering angle for both time intervals after subtraction 
of background. The error bars shown in Figure [2] include 
not only the statistical uncertainty but also the 1.8 % of 
systematical one as estimated the above subsection. 

At the first step, we investigated the polarization de- 
grees, II, and the polarization angle (4> p ), measured 
from north with anticlockwise direction, separately for 
Interval-1 and -2. Then, the response of GAP for irradi- 
ation from 20.0 degree off-axis is modeled by the Monte- 
Carlo method with Geant 4 simulator. The gray solid 
lines are the best-fit functions for Interval-1 and -2. The 
best values are IT = 25 ± 15 % with cf> p i = 159 ± 18 de- 
grees for Interval-1 and II2 = 31±21 % with (f> P 2 — 75±20 
degrees for Interval-2, respectively. Hereafter the quoted 
errors are 1 a confidence for two parameters of interest. 
The significance of polarization detection is rather low of 
95.4 % and 89.0 % for Interval-1 and -2, while the dif- 
ference of polarization angles is significant with 99.9 % 
(3.5 a) level. 

In the next step, we performed a combined fit for the 
two intervals, assuming that the polarization degree for 
Interval-2 is the same as that for Interval-1. This means 
that we treat II as one free parameter to improve the 
statistics with the reduction of model parameters. Here 
the two polarization angles were still free parameters for 
both intervals, because the change of angle is apparent. 
The best-fit polarization degree is n = 27 ± 11 % with 
X 2 = 21.8 for 19 degrees of freedom. We show a A% 2 
map in the (n.,<f> p ) plane in Figure [3l where we represent 
4> P = 4>pi ■ The significance of detection of polarization is 



99.4 % (2.9 a) confidence level. 

Gamma-ray detectors aboard RHESSI and INTE- 
GRAL have complex geometry, the polarization anal- 
ysis may be highly affected by t he instrumental sys- 
tematics as the autho rs mentioned (jKalemci et al.ll2"007t 
iMcGlvnn et al.ll2007t ). In contrast to them, GAP is de- 
veloped for the purpose of GRB polarimetry with the 
high axial symmetry in shape. The systematic un- 
certainty is quantitatively estimated and found to be 
as small as 1.8 % of the total polarization signals, as 
described above. Of course, GAP is well calibrated 
for p olarized gamma-rays before launch (|Yonetoku et al.l 
120 111) . Therefore, for the prompt emission of GRB, we 
conclude that this is probably the most convincing de- 
tection of polarization degree and angles so far. 

4. DISCUSSION 

The emission mechanisms of GRBs which are being 
actively discussed with the spectral and lightcurve in- 
formation collected so far includes synchrotron emission 
of electrons at the optically-thin regions of rclativistic 
jets and quasi-therm al emission from the photospheres 
of jets (for a review. fMeszarosl feOOS). Our reliable po- 
larimetric observations may help further constrain the 
emission mechanisms as well as the structure of the emis- 
sion sites and the magnetic fields. The main results of 
our polarimetric observation of GRB 100826A are: (1) 
The polarization degree averaged over the burst duration 
is II = 27 ± 11 % in the energy range of 70 - 300 keV 
below E p ~ 600 keV; and (2) The PA significantly varies 
from Interval-1 to -2. 

It has been theoretically shown that sizable net polar- 
izations can be obtained both in the synchrotron and 
photospheric emission mechanisms. The synchrotron 
mechanism involves several different types of models 
with respect to magnetic field structure (for reviews, 
lLazzatil[2006t PTbma et a l. 2009); (a) synch rotron in glob- 
ally ordered toroida l field (SO model: jGranotl 120031 : 
ILvutikov et all [2001 : IZhang fc Yaiil l201lh . (b) in ran- 
dom fields on plasma skin d epth scales (SR model; 
lGranotl[2003l:lNakar et aTJl2003[) . and ( c) in random fields 
on hydrodynamic scales (S H model; llnoue et"ail 120111 : 
iGruzinov &: Waxmanl [l999). In the SR model, the mag- 
netic field directions are assumed to be random not 
isotropically but mainly in the plane perpendicular to 
the local expansion direction. In this case a high degree 
of polarization can be obtained only when the observer 
sees the jet slightly off-axis, i.e., 9 V ~ 9j + r — 1 , where 
9 V is the angle between the line of sight and the jet axis, 
and 9j and T are the opening half angle and the bulk 
Lorentz factor of the jet, respectively. On the other hand, 
(d) in the photospheric emission model (Ph model), we 
also expect the strongly polarized radiation when the 
radiation energy is comparable to or sma ller than the 
baryon kinetic energy at the photosphere (|Beloborodovl 
120111 ). Around the photosphere the photon distribution 
is highly anisotropic, which can lead to high levels of 
pol arization through the last electron scatterin gs (see 
also lEichler fc Levinsonll2003l : lLazzati et aLll2004 for the 
other Comptonization models). We need the condition 
8 V ~ Sj + r _1 also in this model to have a high net po- 
larization. The above four models can be consistent with 
our result (1). 

Our result (2) provides us with very important infor- 
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mation, which excludes the axisymmetric jets both in 
the synchrotron and photospheric emission models. If 
the emitting surface is symmetric around the jet axis as 
assumed by most of authors in the SO, SR, and Ph mod- 
els, the direction of the net polarization in the high II 
case is determined to be either parallel or perpendicular 
to the direction from the jet axis towards the line of sight. 
Therefore we conclude that non-axisymmetric structure 
of the brightness and/or the magnetic fields on the ob- 
servable angular scale ~ T" 1 are required to have the 
subst antial change of PA (see also lLazzati fc Begelmanl 
f2009h . 

A possible solution for the results (1) and (2) is the 
patchy emission in the SO, SR, or Ph model, or the 
patchy magnetic field structure, i.e., the SH model. For 
instance, the SH model produces the net polarization as 
II ~ 70% / yN, where N is the number of independent 
patches with coherent magnetic field in the observable 
region. The polarization degree for N ~ 2 — 50 can be 
consistent with the individual results for Intervals- 1 and 
-2 within the 1 a confidence. The angle of the averaged 
polarization is statistically random, so that they can be 
consistent with our result (2). Similar arguments may 
be made in the patchy SO, SR, and Ph models. In the 
patchy SO model, the emission from each patch can have 
different PA easily when 9j ~ T" 1 . In the patchy SR and 
Ph models, we no longer require the off-axis viewing of 
the jet, and the emission from each of off-axis viewed 
patches has high II and different PA. 

How can we further constrain the models with the anal- 
ysis results of other bursts detected by GAP and the 
future polarimetry missions? We note that the statis- 
tical analys i s resu lts of the SO, SR and Ph models in 
iToma et al.l (|2009f) are not applicable since they assume 
the axisymmetric jets. In the Ph model, the patch scales 



are constrained as 9 p > T -1 . Then high local polariza- 
tions are obtained from patches with 9 V ~ 6 p + T^ 1 > 
2r _1 , although their brightness is much smaller than 
that from patches with 9 V < because of the rel- 
ativists beaming effect. Thus it is not expected that 
all the bright bursts have a high net polarization in the 
patchy Ph model. The statistical studies of the GAP 
bursts might provide us some implications for further 
constraining the emission mechanism. 

This polarimetric observation for GRB 100826A has 
been performed in the energy range below E p . In real- 
ity, the spectral indices as of bursts below E p are not 
fully explained in the framework of synchrotron or photo- 
spheric mechanism. In the synchrotron mechanism, the 
radiative cooling of emitting electrons should not be so 
rapid to produce the spectrum sufficiently hard below 
E p , so that the emission region could be thick, unlike 
the thin shell emission assumed in many papers (e.g., 
lAsano fc~T crasawa 2009|). The photospheric models in- 
clude an idea that the emission below E p is superposition 
of many quasi-thermal c omponents with d i fferent tem- 
perat ure and/or V (e.g., IRvde et all 120101 : IToma et"aLl 
l20ll . Clearly, more studies on the energy dependence 
of polarization (as well as the non-axisymmetric struc- 
ture of the jet) in each model are needed to discuss the 
observational results quantitatively. 
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Fig. 1. — Lightcurve of the prompt 7-ray emission of GRB100826A detected by the GAP. We divide the data into Interval-1 and -2 for 
the polarization analysis. 
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Fig. 2. — Number of coincidence 7-ray photons (polarization signals) against the scattering angle of GRB 100826A measured by the 
GAP in 70-300 keV band. Black filled and open squares are the angular distributions of Compton scattered 7-rays of Interval-1 and -2, 
respectively. The gray solid lines are the best-fit models calculated with our Geant4 Monte-Carlo simulations. 
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Fig. 3. — A A.\ 2 map of confidence contours in the (n, <j> p ) plane for GRB 100826A, obtained by the combined fit of the Interval-1 and -2 
data. Here (f> p is the phase angle for Interval-1. The white dot is the best-fit result, and we calculate A\ 2 values relative to this point. A 
color scale bar along the right side of the contour shows levels of A\ 2 values. The null hypothesis (zero polarization degree) can be ruled 
out with 99.4 % (2.9 a) confidence level. 



